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Abstract 
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In this paper we have solved the hydrodynamic equations governing the spher- 
ically symmetric isothermal accretion (wind) onto (away from) compact ob- 
jects using various pseudo-Schwarzschild potentials.These solutions are es- 
sentially one parameter solutions in a sense that all relevant dynamical as well 
as thermodynamic quantities for such a flow could be obtained (with the as- 
sumption of a one-temperature fluid) if only one flow parameter (temperature 
of the flow T) is given. Also we have investigated the transonic behaviour 
of such a flow and showed that for a given T, transitions from subsonic to 
the supersonic branch of accretion (wind) takes place at different locations de- 
pending on the potentials used to study the flow and we have identified these 
transition zones for flows in various such potentials. 
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1 Introduction 



A number of modified Newtonian potentials of various forms are available in 
literature which accurately approximate some general relativistic effects im- 
portant to study accretion discs around a Schwarzschild black hole. Such po- 
tentials may be called 'pseudo-Schwarzschild' potentials because they nicely 
mimic the space-time around non-rotating/ slowly rotating compact astrophys- 
ical bodies. Recently we have shown that^ (Das c& Sarkar, 2001, hereafter 
PI) though the potentials discussed in this paper were originally proposed to 
mimic the relativistic effects manifested in disc accretion, it is quite reason- 

^Send any correspondense/ offprint requests to Tapas K. Das {tapas@iucaa.ernet.in) 
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able to use most of these potentials in studying various dynamical as well as 
thermodynamic quantities also for spherical Bondi-type^ accretion of adiabatic 
fluid onto Schwarzschild black holes. In PI, the space variation of various dy- 
namical and thermodynamic quantities had been studied for a two parameter 
adiabatic accretion and wind system i.e., flow parametrized by the specific en- 
ergy and the polytropic constant of the flow. In this paper, we would like to 
continue our investigation of spherical accretion for the same set of pseudo 
potentials used in PI but for isothermal accretion (wind) which may be charac- 
terized by a single parameter, i.e., the flow temperature which can be approx- 
imated by the proton temperature for an one-temperature fluid system. 
Owing to the fact that, close to a spherically accreting Schwarzscild black hole, 
the ratio of electron number density to the photon number density is propor- 
tional to ^/r (where r is the radial distance measured from the central accretor 
in units of Schwarzschild radius Vg = "^^^.f" ), a lesser number of electrons 
would available per photon with decrease of r, and close to the accretor, mo- 
mentum transfer by photons on the accreting fluid (or wind) might be an effi- 
cient process which may keep the flow temperature roughly constant (at least 
upto the sonic point) and hence isothermality may not be an entirely unjusti- 
fied assumption^. However, this efficiency for momentum transfer falls of for 
large r and far away from the black hole isothermality assumption may brake 
down. 



2 Governing equations 

We take the following pseudo potentials: 
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to study the flow. Whereas and $2 had been proposed by Paczynski and 
Wiita^ (1980) and Nowak and Wagoner^ (1991) respectively, $3 and $4 are 
suggested by Artemova et. al.^ (1996) (see [6] and PI for detail discussion 
about these potentials). Hereafter, we will denote any ith potential as $i where 
{i = 1, 2, 3, 4} corresponds to {$1, $2, ^3, ^4} respectively. 
For isothermal equation of state. 
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where P, p, T and Cg are the pressure, density, constant temperature and the 
constant isothermal sound speed of the flow respectively; R and /i being uni- 
versal gas constant and the mean molecular weight. Cg and T can be related 
as. 



a = @T 



(3) 



2 



Where 9 = a constant, uih ~ rup being the mass of the Hydrogen atom. 
Integrating tixe radial momentum conservation equation 

du du 1 dP ^ / ^ 
dt dr p dr * 

and the equation of continuity 

using eqn. (2) and (3), we obtain two conservation equations for our flow as 
(in system of units used in PI): 



u: 



+ STlnpi + $i = C 
M = A-KpiUir^ 
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Where C is a constant, u and M being the dynamical flow velocity and the 
mass accretion rate respectively. Subscript i indicates that respective quanti- 
ties are measured for a particular The space rate of change of dynamical 
velocity can be written as: 
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dr 
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where ^[ = From eqn. (4b), one can easily obtain the sonic point condi- 
tions as. 
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where <I> , 



is the value of ^ at the corresponding sonic point r^. ^ at the 



sonic point may be obtained by solving the following equation 
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where 



is the value of at sonic point. Mach number of the flow Mj 
for a given T and can be expressed as. 

Mi = Ui@-^-^T-^-^ = 



iTMpi-'^r-'^ + 2$iG-iT-i + (2T2 - 1) Inpi 
where M is the mass accretion rate defined in eqn. (4b). 



(5) 
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3 Results 



For a given flow temperature T, one can solve eqn. (4d) to find out the sonic 
point rl for any particular ith potential In figure 1 we show the varia- 
tion of sonic points (for #j=i_4) with constant flow temperature. We observe 
that for all $j=i_4's, sonic point anticorrelates with T non-linearly and mono- 
tonically. This is obvious because as sound speed Cs is constant throughout 
isothermal flow, it is only the dynamical velocity u whose change would con- 
trol the change of Mach number, and as Cg is proportional to VT, the more 
will be the flow temperature, the more will be Cs and the less will be the Mach 
number M (M = ^) for a particular radial distance and the less will be the 
location of rc where u becomes equal to Cs, i.e., M becomes unity. It is evident 
from the figure that for a given flow temperature (fixed T): 

rl > rl > rl > rl (6) 

For a fixed T, it is now quite straightforward to simultaneously solve eqn. (4a- 
4e) and eqn. 5 to get the integral curves of motion for all <I>i=i_4's to study the 
transonicity of the flow. In figure 2, for T = 5 x lO^'^ °K, we plot all the integral 
curves of motion for accretion as well as for all $j=i_4's described in §1. For a 
given flow temperature and for a particular distance r = Tq, if we define to 
be the Mach number attained at Tq for a specific <I>i, it is clear from the figure 
that 

Ml > > Ml > Ml (7) 

while the sequence is just reversed for the wind. As velocity profile for the 
isothermal flow is identical in nature of its Mach number profile with a scale 
shift of (6r)~°'^, one can conclude that produces the steepest velocity gra- 
dient for accretion and flattest velocity gradient for wind. Combining eqn. (6) 
with (7), one can interpret that for isothermal accretion, Paczyhski and Wiita"^ 
(1980) potential $i would produce the maximum spatial rate of change of ki- 
netic energy of the accretion for a given T while the same is produced by po- 
tential proposed by Nowak and Wagoner^ (1991), i.e., $2/ for wind. Similarly, if 
for any one defines the spatial gradient of Mach number to be which 
might be regarded as the measure of the 'transonicity' of the flow (because 
more will be the value of the faster the flow will become supersonic, i.e., 
the further will be the sonic point away from the black hole), we find that for 
accretion branch. 



-g^ , e.g., the absolute value of increases non-linearly 
and monotonically as the flow moves towards the black hole and sequence for 



dr 
dM^ 



follows the same trend shown in eqn. (7). However, for wind branch, 
does not monotonically increase as the flow moves away from the accre- 
attains a maximum value in the subsonic fart of the wind and 



dr 

tor rather -g^ 

then starts falling non-linearly, see figure 3. The location of the maxima can be 
obtained by solving the following non-trivial equation for (the location of 
the peak): 

k)' *;,(rp, Up) + rl^^irp, Up) - 2eT = (8) 
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Where *'(r, u) and ri'(r, u) are two complicated functions of r and u defined 
as: 



The superscript i indicates the value obtained for a particular ith potential. 
Also we obtain that the location of the peak maintains the following sequence: 



Thus we successfully analyze various aspects of transonic behaviour of the 
spherically symmetric, isothermal accretion characterized by its flow temper- 
ature T. 

Acknowledgements 

We are thankful to Prof. I. D. Novikov and Prof. P. J. Wiita for useful dis- 
cussions. 

References 

[1] T. K. Das and A. Sarkar , Astron. Astrophys. 374, 1150 (2001) (PI) 

[2] H. BondiMonthly Notices Roy. Astron. Soc. 112, 195 (1952) 

[3] T. K. Das, Monthly Notices Roy. Astron. Soc. 318, 294 (2000) 

[4] B. Paczyhski and P. J. Wiita, Astron. Astrophys. 88, 23 (1982) 

[5] A. M. Nowak and R. V. Wagoner, Astrophys. J. 378, 656 (1991) 

[6] I. V. Artemova, G. Bjomsson and I. D. Novikov Astrophys. J. 461, 565 (1996) 





' p ' p ' p ' p 



(9) 



5 



Figure Captions 



1. Fig. 1: Variation of sonic point with flow temperature for various $i=i_4's. 
While flow temperature T (in units of Tio = 10^*^ °K) is plotted (in loga- 
rithmic scale) along x axis, sonic point (in logarithmic scale) is plotted 
along y axis. 

2. Fig. 2: The integral curves of motion drawn for various $i=i-4. While 
the radial distance from the black hole (in units of Vg) is plotted along x 
axis in logarithmic scale, Mach number M of the flow is plotted along y 
axis. The sonic points for various $j=i_4's are obtained as {r^ , , , } = 
{28.05, 23.14, 26.6, 27.09} Vg respectively. 

3. Fig. 3: Variation of gradient of Mach number for a given T with respect 
to spatial distance. While the radial distance from the accretor (in units of 
rg) is plotted along x axis in logarithmic scale, the spatial rate of change 
of Mach number (^^^ is plotted along y axis. It is clear that the curves 
obtained for all $i=i-4 produces a 'peak' location of which is designated 
as Vp in eqn. (8). is different for different ^/s for a particular given 
flow temperature T. 
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Fig. 3: 
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